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1—3 ■ Abstract 

New structure function measurements from fixed target experi- 
ments and especially HERA are reviewed. The extraction of parton 
I ' distributions from these measurements is discussed with special em- 

phasis on systematic problems. New information from Drell-Yan and 
. direct photon production experiments are also presented . Finally the 

| present uncertainties of our knowledge on parton distributions and on 

a s from DIS experiments are discussed []. 
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1 Introduction 

V ■ 

Ph! Our knowledge of nucleon structure functions and of the parton distributions 

(PDF's) derived from them has steadily improved due to both the improve- 
>• . ment of a large variety of measurements and a more sophisticated theoretical 

^ ■ treatment of hard scattering processes in perturbative QCD. Structure func- 

c3 ■ tions and PDF's are needed for two reasons: 

• They are a necessary input for all hard scattering processes involv- 
ing nucleons to make precise predictions in the standard model and 
of course to look for deviations from theses predictions, and to make 
predictions for signals and backgrounds at colliders especially the LHC. 

• They contain important information about the underlying physics of 
hadrons and allow stringent tests of perturbative QCD. 

A hard scattering process in a hadron hadron collision is shown in figure 
|l]. Two partons with momentum fractions X\ and x<i scatter and produce a 
1 Invited talk XVIII Physics in Collision 



1 



Figure 1: Example of a LO hard scattering process in hadron-hadron scatter- 
ing. 



heavy state with mass M which can be jets, 7+jet, heavy Bosons or Quarks 
or new particles like Higgs or SUSY. The differential parton-parton cross 
section da^ is given by pQCD, the total cross section is given by 



where the sum goes over all parton flavours which contribute. In LO pQCD 
the cross section factorises into the parton luminosity which depends on 
the parton densitites /j, fj and the differential parton-parton cross section. 
The parton densities are therefore universal to all hard scattering processes 
provided they are extracted in higher order (so far in NLO) and corrected to 
the leading order diagrams. The parton densities depend on the fractional 
momentum x and the typical scale fi 2 of the process. 

Deep inelastic lepton nucleon scattering (DIS) is the basis for our knowl- 
egde of PDF's. The kinematic plane of available measurements in (x,Q 2 ) 
from fixed target and HERA experiments is shown in figure 0. 

The data cover now a huge region 



where HERA has added 2 orders of magnitude in both variables. The kine- 
matic plane for hard scattering processes at LHC (14 TeV) is compared in 
figure | where the mass M of the produced system is used as scale. Also 
shown are lines of constant rapidity y of the system M which is proportional 
to (xi — X2- It is obvious that due to the high CMS energy of LHC a very 
large fraction of interesting LHC physics is physics at low x ( here taken as 
x < 10~ 2 ) JJj. For example the production of Higgs particles in the mass 



da = ^2 dx 1 dx 2 f i (x 1 , ^ 2 )fj(x2, /i 2 ) * do- ij (p 1 ,p 2 ,a s (n 2 ), M 2 / fi 2 ) 




2 * 10" 5 < x < .75 ; 1 GeV 2 < Q 2 < 40000 GeV 2 
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Figure 2: Available DIS data from fixed target and HERA. 



range 100 GeV < Mh < 500 GeV in the rapidity range \y\ < 2.5 requires 
the knowledge of the gluon density for x- values as low as 10~ 4 which is now 
available due to HERA. 

In order to apply our present knowledge on PDF's we have however to 
extrapolate the parton distribution by up to 3 orders in magnitude in Q 2 . 
This is possible using the pQCD evolution equations. At large x the DGLAP 
equations are expected to provide a sufficiently good approximation. This is 
however rather doubtful at small x because terms proportional to a s *lnl/x 
have been neglected in the derivation of DGLAP compared to the a s * InQ 2 
terms which have been summed. This does not look justified at small x. One 
important question at HERA is therefore the test of QCD dynamics at small 
x. 
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LHC parton kinematics 
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x 12 = (M/14TeV) exp(±y) 
Q= M 
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Figure 3: Kinematic plane of hard scattering processes at LHC. Lines of 
constant rapidity y for the produced system with mass M are given as dashed 
lines. 
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2 New DIS data 
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Figure 4: The structure function vs. Q 2 for fixed values of x for recent 
HERA and for fixed target data. The solid lines are the result of a NLO 
pQCD fit to these data by HI as discussed below. 



Progress in the measurement of nucleon structure functions is steady but 
rather slow. It takes long time and hard work to make reliable measurements 
and to study and reduce sytematic effects. Nevertheless the best measure- 
ments are all dominated by correlated systematic errors which makes it diffi- 
cult to analyse these data in a consistent quantitative way. From fixed target 
experiments two important new results have been added. The NMC collabo- 
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ration has published the final analysis of their data taken back in '89 with a 
small angle trigger |2| . These data enlarge the low x range and therefore the 
overlap with HERA and also allow to measure R = (Tl/^t- The CCFR neu- 
trino collaboration at Fermilab has published an improved analysis of their 
data first published in '93 ||. They have essentially improved the energy 
calibration for hadrons and muons. Their structure function measurements 
have a significant impact on our knowledge of a s , to be discussed at the end. 
The HI collaboration has made a very substantial step towards large x and 
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Figure 5: Fixed target and HERA measurements of F 2 P for a fixed value of 

Q2 = 15 Q eV 2 



1.2 



CCFR • 


O.Q075 


NMC ■ 


O.D08 


E665 a 


0.00893 i 


H1 T 


0.008 


III' 


HIprelA 


0.008 1 n i 

1 ttl 
.if" 




* 













1 10 io 2 

Q ! (GeV") 

Figure 6: Fixed target and HERA measurements of for a fixed value of 
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Q 2 which allows for the first time to directly test the Q 2 -evolution at large 
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x over two orders of magnitude with reasonable precision. Also more precise 
data in the low x region has been provided by both HI and ZEUS. A compila- 
tion of some of the most precise measurements of the structure function F 2 ep 
is shown in figure |] versus Q 2 for fixed values of x. Figures |5| and ^ show all 
available data points for a fixed value of Q 2 = 15 GeV 2 resp. a fixed value 
of x=.008 in order to see the relative precision and the matching of fixed 
target and HERA data. Figure [5] demonstrates very clearly one of the most 
important contributions of HERA so far: The structure function at small x 
and therefore the qq sea rises by more than a factor 3 in the HERA range 
towards small x. This has severely affected the predictions for processes at 
LHC both for pp and heavy ion collisions. The structure function for fixed 
x=.008 rises by a factor 5 with Q 2 demonstrating huge scaling violations 
whch have to be explained by pQCD. The overall agreement among the two 
HERA experiments and the matching with fixed target experiments is good. 
Our knowledge on F% p has presently a systematic uncertainty of about ±3% 
over the whole x-range 10~ 4 < x < .65 mainly limited by the uncertainties 
of energy scales. 

3 Analysis of the low x region and the gluon 
distribution 

The low x data on F 2 show large scaling violations . For large enough Q 2 
these should be described by the pQCD evolution equations: 

dZ(x,Q 2 )/dlnQ 2 = a s /2n (P gq ® g + P qg <g> E) 

dxg{x, Q 2 )/dlnQ 2 = a s /2-K (P gg ®g + P qg <g> S) 

Here £ = J2f x (lf J rQf) is the sum of all quark and antiquark momentum dis- 
tributions. The parton distributions evolve with Q 2 because quarks radiate 
gluons, gluons split into quar- antiquark pairs and gluons split into gluons. 
At small x the gluon radiation from quarks can be neglected. The evolution 
is then dominated by the evolution of the gluon distribution. The splitting 
functions Pij depend on terms proportional to 

a s -k InQ 2 , a s * lnQ 2 lnl/x and a s * Inl/x 
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Figure 7: Measured slopes dF 2 /dlnQ 2 obtained from linear fits in InQ 2 to the 
HI data. The solid line gives the slopes as obtained from the NLO Q CD fit. 



If terms proportional lnl/x are neglected, then we arrive at the DGLAP 
evolution equation, if instead terms proportional to lnl/x are kept instead of 
InQ 2 terms we arrive at the BFKL equation. 

The first thing to do is therefore to check the evolution directly. Figure 
[7| compares the measured slopes dF 2 /dlnQ 2 (x) (by linear fits in InQ 2 for 
fixed x ) with the predictions of a DGLAP fit to the data ( HI fit discussed 
below). There is good agreement down to x — 10~ 4 . This illustrates, to our 
big surprise, that the DGLAP aproximations works even at very small values 
of x. Other attempts to see 'BFKL' like effects in the final state , which are 
not discussed here, were also negative. 

This observation gives us the possibility to measure the gluon distribution 
from the measured slopes dF 2 /dlnQ 2 (x). To rather good approximation the 
relation at small x is given by: 

dF 2 /dlnQ 2 (x) « 10/27vr * a s (Q 2 ) * xg(2x, Q 2 ) 

The HERA experiments therefore measure the gluon distribution at small x 
which is of prime importance for physics at LHC. 

Both HERA collaborations have determined the gluon distribution from 
a NLO DGLAP fit to their data (Q 2 > 1.5 GeV 2 ) in combination with ep and 
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HI (Ql - 1 GeV 2 ): 



x 9( x >Qo) = A g x B a(l — x) c 9 «— 3 Parameters 

xUv(x,Ql) - A a x B "(l - a;) c "(l + D u x + E u ^) 

xdv{ X ,Qfo = A d v B *{l - x) c *{l + D d x + E d *Jx~) 

xS(x,Q%) = A s * Bs (l -m) Cs (l+ D s x+ Egvfa) 

m assume u = d 

• stiange quark 20% of total sea 

• a s (M z ) - 0.118 

Figure 8: Parametrisations of parton momentum distributions as used by HI 
and choice of parameters. The solid line gives the slope as obtained from the 
NLO QCD fit. 



ed data from BCDMS and NMC. This is explained below for the example 
of HI || which uses 15 free parameters plus some relative normalisation 
parameters. The resulting gluon distribution is shown in figure |9] at Q 2 = 
20 GeV 2 together with the corresponding determination from ZEUS Q. The 
gluon density rises strongly towards small x and xg(x) reaches values of 25 at 
x = 10 -4 , a factor 15 larger than the sea quark density. This underlines the 
predominace of the gluon evolution by gluon splitting whereas the splitting 
of a gluon into a qq pair is relatively rare. The error bands given include 
all known strongly correlated systematic errors. The gluon distribution in 
the x- range between 10~ 2 and 10~ 4 is therefore known to about 15% from 
HERA. It should be noted that the charmed quark contribution to F2 is very 
substantial (up to 20%) at small x. This contribution is calculated directly 
in NLO pQCD || via the photon gluon fusion process. This contribution has 
been directly measured by both HERA experiments using DIS events with 



D*'s. The preliminary result of ZEUS is shown in figure |10| together wih the 
pQCD prediction based on the gluon distribution of figure |9|. There is good 
agreement showing the consistency of the procedure. 
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Figure 9: Resulting gluon distribution for Qq = 20GeV 2 from HI, ZEUS and 
NMC NLO pQCD fits. The error bands include the estimate of correlated 
systematic errors. 
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Figure 10: F^ iarm as measured by the ZEUS collaboration using D* produc- 
tion in DIS compared to the NLO prediction based on the 7 — gluon fusion 
process. 
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4 Transition to low Q 2 at HERA 



So far we have discussed the 'DlS-regime' where we beleive pQCD is appli- 
cable without special justification. It's worthwhile to study the transition 
region from the DIS regime to the region down to Q 2 = , e.g. 'photo- 
production'. This has been possible at HERA due to the installation of 
special low angle electron calorimeters which allow to measure inelastic neu- 
tral current scattering down to Q 2 = 0.11 GeV 2 M. In addition we can 



• ZEUS BPC 1995 ■ ZEUS SVTX 95 Prel. • E665 

- DL t HI SVTX 95 • ZEUS, HI yp 

■ GRV(94) a ZEUS 94 
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Figure 11: The total 7*p cross section vs. W 2 for low values of Q 2 compared 
to DGLAP (GRV) and soft pomeron regge (DL) predictions. 



of course measure the total photoproduction cross section aq?{W 2 ,Q 2 = 0). 
The theoretical concepts used to describe DIS and photoproduction can be 
easily related. The structure function F 2 (x, Q 2 ) and the 7p cross section are 
related by: 

a^ p (W 2 ,Q 2 ) ~ 4na 2 /Q 2 F* p (x,Q 2 ) 

W 2 p ~ Q 2 /x at low x 

Scattering of electrons at low x is therefore equivalent to high energy 7^ 
scattering. Since photons behave like hadrons if they interact with nucleons 
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it should not be surprising that the concepts which are used to describe high 
energy hadron-hadron scattering are suitable to discuss the low x region. 
Figure shows crj p as funtion of W 2 for different values of Q 2 . Its well 
known that the photoproduction cross section has the same rise with W 2 
as observed in p-p scattering: aj?(W 2 , Q 2 = 0) ~ (W 2 ) x where A = 0.08 is 
the intercept of the 'soft Pomeron trajectory'. The HERA data at low Q 2 
show that the increase of o~t with W 2 becomes larger with increasing Q 2 , 
therefore A rises with Q 2 as shown in figure [E| from a value close to .08 at 
low Q 2 to A ~ .35 at Q 2 = 40GeV 2 . The measurements are compared to 

fa 0.5 r 
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Figure 12: The exponent A /rom a fit F 2 ~ x A to t/je i/i dato at small x vs. 
Q 2 - 



a soft pomeron model of Donnachie and Landshoff (DL, dashed curves) and 
a pQCD fit based on the DGLAP evolution (dash-dotted line). We see a 
smooth transition from regge models to the DIS description around Q 2 = 
.8 GeV 2 . Figure p~3| shows another interesting aspect of the transition region. 
The slopes dF2/dlnQ 2 have been directly determined from HERA data as a 
function of x. These slopes rise towards low x from x = 10 _1 to x = 10 -4 
in agreement with DGLAP predictions given by the triangles but then fall 
again towards even lower x. The slopes at lowest x are in agreement with 
regge fits based on the 'soft Pomeron' concept. The transition occurs around 
x = 10" 4 at Q 2 sa lGeV 2 . 
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Figure 13: Slopes dF 2 / 'dlnQ 2 measured in the transition region compared to 
predictions of a DGLAP fit and of a Regge fit. 
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5 Large x and large Q 2 



This is the kinematic region where we want to test the standard model by 
looking for e.g. contact interactions or new massive particle production (see 
talk by J. Meyer at this conference). The questions here are; 

• how reliable is the standard model prediction e.g. how well can we 
extrapolate parton densities from fixed target experiments to HERA. 

• what are the experimental limitations for a precise measurement at 
HERA. 



The preliminary high x data from HI is shown in figure [L4| where the reduced 
cross section o ~ F 2 up to electroweak interference effects is shown versus 
Q 2 for different bins in x compared to BCDMS and SLAC data. Also shown 
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Figure 14: The reduced cross section for large values of x vs. Q 2 compared to 
two DGLAP fits to data with Q 2 < 120 GeV 2 and inclusing the large x data. 



is a QCD fit which used only low Q 2 data. The HI measurements at an 
average < Q 2 >= 10000 GeV 2 have an uncertainty of about 14% limited 
mostly by the uncertainties in the energy calibration. We therefore have so 
far directly checked the Q 2 evolution over 2 orders of magnitude to that pre- 
cision. Further improvements will require better calibration in addition to 
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higher statistics. This is indeed possible at HERA. The energy calibration 
can be directly obtained from the data by using the redundancy in the kine- 
matic reconstruction. At HERA we measure angle G e and energy E e of the 
scattered electron as well as angle Qh and energy Eh of the hadronic system 
(the quark) . The energy of the electron can therefore also be calculated from 
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Figure 15: Comparison of the HI calibration of the electron energy in the 
liquid argon calorimeter compared to the energy form the double angle method 
vs. the calorimeter position along the beam axis. 



the two angles alone Eda = /(@e, The energy scale of the electromag- 
netic calorimeter vs. polar and azimuthal angle was therefore determined by 
comparison to the double angle energy. This comparison is shown in figure 



15. This procedure determines the energy scale for the scattered electron to 



about 1% for O e > 15° and to about 3% for smaller angles where the 3% are 
limited by statistics. 

HERA will substantially increase its luminosity in 2000. High statistics 
data will then allow better calibration and therefore very good and reliable 
checks for deviations of the standard model. If a deviation should be ob- 
served due to new physics then we are pretty sure that such an effect could 
be established beyond doubt in contrast to similar measurements at hadron 
colliders because only quarks can contribute and because the energy calibra- 
tion is provided by the data themselves. 
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6 Global fits to parton distributions 



The determination of a complete set of parton distributions from available 
data is the domain of two groups of phenomenologists, the CTEQ group 
centered at Fermilab and the MRS group at Durham. They provide libraries 
of particle density functions (PDF) which can be easily used and which are 
regularity updated if new data become available. 

Their global analyses use fit methods very similar to the one described 
above for HI. They use however additional input in order to also separate 
the flavour content of valence and sea quarks and to better constrain the 
gluon distribution. Main additional input comes from the following sources: 

• Neutrino measurements from the CCFR collaboration || which deter- 
mine 

-the shape and magnitude of the valence distribution x(u v + d v ) from 
the structure function xF 3 

- the magnitude of the strange sea s(x)/(u(x) + d(x)) ~ 0.25 from sin- 
gle charmproduction v + s — > fi + c. - a s determination from scaling 
violations of the nonsinglet structure function. 

• Data to determine u(x)/d(x) mostly in the valence region 

- The direct NMC measurement of u(x)/d(x) and the measurement 
of the asymmetry in W-production by CDF which is the most sensitive 
measurement at large x M. 

• The determination of u/d by Drell-Yan experiments 

• Experiments measuring direct hard photon production in pN scattering 
as a direct measurement of the gluon distribution at large x. 

No attempt is made here to go into details. Comprehensive summaries can 
be found in references || and HlQfl . The results which are summarised here 



in order to illustrate sucesses and problems of this approach are taken from 
the recent fit of the MRS group using the PDF set MRS98 0. This basic 
PDF set fixes the value of a s (mz) = -118 to the world average. 

Global fits have the problem that data sets are incompatible with each 
other and that practically all data sets are systematically limited with highly 
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correlated errors. There is no chance to include these errors - even if they 
would be fully available - into the fits. This may lead to systematic biases 
of the results and it makes it practically impossible to evaluate reliably the 
uncertainties of the parton densities. 



6.1 DIS data in the global fits 




Figure 16: Measurement of F% N (x,Q 2 ) by the CCFR collaboration compared 
to MRS98 prediction. Only measurements with x > 0.1 have been used in the 
fit. 



There is overall reasonable agreement. Some data sets give rather large 
contributions to x 2 I dof like e.g. the BCDMS ep and ed and the ZEUS 
structure functions but the fits disregard the systematic errors. Both MRS 
and CTEQ have decided to use single charm production from the CCFR 
neutrino experiment to fix the strange sea. This however has a problem as 
illustrated in figure [16|. The measured structure function F% N in the sea 
region ( x < 0.1) is 10 % to 20 % higher than expected whereas the NMC 
muon data is well described. Obviously one has to make a choice and MRS 
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has decided NOT to use the low x CCFR data in the fit. There is however 
no explanation for this discrepancy. 

6.2 Drell-Yan processes and the u/d ratio 



E772 (p d -f^VX) PLAB = 800GeV 
10 ». 




0.08 0.09 01 -J x 0.2 03 0.4 



Figure 17: Drell-Yan cross section of experiment E772 for the process 
pd — ► fi + fi~X j71|/ compared to the NLO predictions based on the MRS98 
parametrisation 

Lepton pair production in hadron-hadron scattering outside the \l/ and Y 
mass ranges as well as the production of the vector bosons at the pp collider 
are rather well described by the Drell-Yan process. Figure |I7] shows the 
differential cross section for the reaction pd — > fi + fi~X as measured by the 
E772 experiment at Fermilab at 800 GeV laboratory energy [JTTJ compared 
to NLO pQCD predictions based on the MRS98 parametrisation. 

The overall agreement is quite good, at a level of about 10%, except in 
the region of large xp and small \Jr where systematic deviations of up to a 
factor 2 are observed. This region corresponds to scattering processes where 
the first parton carries a momentum x\ ~ 0.4-=-0.6 whereas the second parton 
is at small x<± ~ .02. Both parton densities should therefore be well known 
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Figure 18: The cross section ratio a(pd — > fi + n~X) /2 * o~(pp — > fi + fi^.X) as 
measured by experiment E866 fiTBrf compared to predictions based onu/d = 1. 
and on the MRS98 PDF with d> u 



such that the origin of the discrepancy is either an experimental problem or 
a problem in the theoretical calculation. 

If one disregards this kind of problems then Drell-Yan data can be used 
to measure the flavour content of the sea quarks especially the u/d ratio. We 
know for quite some time that there are more down quarks in the sea than 
up quarks because the Gottfried sum rule as measured by NMC @ 

/ (Fg p - F$ n )dxlx = 1/3 - 2/3 f\d - u)dx = .235 ± .026 
Jo Jo 

differs significantly from 1/3 and because the Drell-Yan experiment NA51 
measured d/u > 1 at x ~ .5 [0. A new dedicated experiment (E866) at 
Fermilab has made very precise measurements in '97 to clarify this question 
flIB| . E866 measured the cross section ratio a{pd — > fi + fi~X)/2 * o~(pp — ► 
fi + fi~X) directly which depends on d/u. This cross section ratio is shown in 
figure [18| versus x compared to predictions for d/u = 1. The MRS98 fit uses 
this data to fix the flavour ratio of the sea quarks. 

The measured difference between down and up quarks in the sea is well 
compatible with models where the nucleon has a pion cloud. 

The ratio of up and down quark densities at large x e.g. in the valence 
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region can be derived from F% n jF% as measured by NMC 0. A very sen- 
sitive new measurement at high scale has been provided by CDF || which 
measured the asymmetry of the decay leptons from bosons at the pp 
collider. This measurement has been selected by MRS to determine u/d in 
the valence region together of course with the BCDMS structure function 
measurements on proton and deuterium. 



6.3 Measurement of the gluon distribution 

The gluon distribution is notoriously difficult to measure over the whole x 
range. In the small x region (x < .1) the gluon distribution is reasonably well 
measured from dF 2 /dlnQ 2 by the HERA experiments and NMC as described 
above. This does not help at larger x because there the scaling violations due 
to quarks radiating gluons dominate over the gluon splitting contribution. A 
very strong constraint for the gluon distribution at medium and large x is 
however given by the energy momentum sum rule. 

A direct measurement of the gluon distribution at medium and large x 
is very desirable. The best suited process is direct hard photon production 
in pN scattering. This process in leading order is absolutely dominated by 
quark-gluon scattering, qq processes give a relatively small and well known 
contribution. Very precise new measurements for this process over a large 
kinematic range have been provided by experiment E706 as shown in 



figure [Bl The invariant cross section for the process pBe — > 7X at 530 
GeV and 800 GeV has been measured for transverse momenta of the pho- 
ton between 3.5 and 11 GeV/c which covers a range of the gluon fractional 
momentum .3 < x 1 < .5. The NLO pQCD prediction based on the 'canoni- 
cal' gluon distribution is unable to describe this data, it is about a factor 2 
too low. The pQCD predictions have a rather large scale dependence; what 
is however worst is the fact that single photon production requires a large 
'intrinsic' kx ~ 1 GeV/c of the incoming gluon as e.g. shown by the im- 
balance of photon and jet momenta of 7 — jet events. Adding a value of 
k T = 1 GeV I c changes the cross section by about a factor 2 independent of 
Pt H3- Both observations show that the theoretical understanding of this 
process is not sufficient and as long as this does not change this process can- 
not be used to rigoriusly constrain the gluon distribution at large x. This is 
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Figure 19: JTie single photon cross section as measured by experiment E701 
in pBe interactions at 530 GeV. The different curves give NLO predictions 
based on MRS98 using two scales and adding an intrinsic kx = 1.2 GeV . 



rather unfortunate because these measurements can in principle decide be- 
tween conventional and 'unconventional'shapes of the gluon distribution at 
large x which have been proposed to explain the excess of dijet events by 
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6.4 How well do we know the parton distributions? 

The parton momentum distributions from the MRS98 fit are shown in figure 
20| for a scale Q 2 = 20 GeV 2 . These can be used to predict hard scattering 
processes. The question of course is: what are the errors for such a prediction. 
This information is not really available. The procedure to estimate this error 
by comparing different parton density parametrisations which can be found 
rather often is completely inadequate. It happens too often that the addition 
of new data changes the PDF data sets far outside the error estimates given 



before. Examples of recent changes are given in figure ^1] taken from the last 
MRS global fit paper M. 
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Figure 20: Parton momentum distributions for a scale Q 2 = 20GeV 2 for the 
MRS98 fit /I/. 



It shows the relative changes for valence quarks and gluons for the MRS 
global fits of '97 and '98. Most surprising is the change of the down quark 
density at large x by up to 25 %. This is the result of including the CDF 
asymmetry measurements in the fit. 

The knowledge of parton distributions is experimentally limited by the 
fact that practically all data sets are dominated by systematic errors where 
calibration errors are dominant. Different data sets also show inconsistencies, 
such that we have to make a choice which ones to use. New data alone will 
not be sufficient to make progress. We need a better treatment of systematic 
errors in the fits where possible and we have to be more selective in choosing 
the data sets entering the fits. Those parts of data sets which have large 
correlated systematic errors should be disregarded in future. 

On the theoretical side more studies are needed to see if the parametrisa- 
tions have enough flexibility. NNLO calculations for DIS would be welcome 
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Figure 21: Relative changes of partem densities vs. x for MRS fits from ? 97 
and 7 98. 

in order to reduce the scale dependences. Most of all we need however a 
satisfactory theoretical description of single photon production to nail down 
the gluon distribution at large x. 

Our knowledge of parton densities to my personal judgement can be 
roughly summarised as follows 

• the cross sections for qq and qq processes can be predicted with an 
error Aa/a < 10% for 10~ 4 < x < .3 . The error increases to about 
20 % for larger values of x. 

• for gluon-gluon scattering we have Aa/a < 20% for x < .1. The error 
increases to 30% up to x=.3 and is about 60% for .3 < x < A. These 
error estimates are more conservative than what is normally given by 
the authors of PDF's. 

Experimental improvements can be expected from HERA in future for quarks 
and gluon densities at small x with better statistics and improved systematic 
understanding of the detectors. 
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6.5 Determination of a. from DIS. 



The uncertainty of a s gives a non negligible contribution to our knowlege of 
parton densities at high scales. The strong coupling constant a s can be well 
determined from the observed scaling violations in DIS because the pQCD 
prediction does not suffer from large theoretical uncertainties and because 
experimental errors are also small. Since several years the world average of 
a s from DIS experiments was always quoted to be low compared to LEP 
measurements. It should be noted that this low value of a s was enforced 
by only two experiments: the BCDMS muon experiment which published 

|TE| and by the CCFR neutrino 



.113 ± .003 exp ± .004 th . 



a value a s (mz) 
experiment which in '93 published a value a 



\ m z, 



.111 ± .002 ± .003 



syst 



17| The CCFR collaboration has recently published a reanalysis of the same 
data ||. After significant improvements of the calibration for hadrons and 
muons the old published value is changed to a s {m z ) = .119±.002 exp ±.004 t / l .. 
This leaves BCDMS as the sole data set which requires a small value. Its 
worth while to have a closer look to this data. The large x data for F| d from 
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neighbouring data points from SLAC and the MRS98 fit with a s = .118. 



24 



BCDMS is shown in figure ^ together with SLAC data and a global QCD fit 
with a s = .118. For each x bin the data points at low Q 2 show large and very 
significant deviations from the fit curve but also from the SLAC data points 
leading to a very poor x 2 /dof =248/174. The data points with systematic 
deviations are at low y where every DIS experiment has the most serious 
systematic problems due to energy calibration. This is also clearly stated in 



the BCDMS paper |1| which says: .." the data agree with SLAC within the 
large correlated systematic errors..'. A fit which uses these statistically very 
precise data points disregarding the systematic errors will therefore severely 
bias the result and necessarily lead to wrong results. 

My conclusion is that i) DIS data is in agreement with a world average of 
a s ~ .118 and ii) people which make global fits for PDF determination should 
care more about systematics, best in collaboration with experimentalists. 
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